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ABSTRACT 

We consider an accretion-disc origin for the broad and luminous forbidden-line emis- 
sion observed in ultraluminous X-ray (ULX) sources CXOJ03383 1.8-352604 and XMMU 
122939.7H-075333 in globular clusters hosted by elliptical galaxies NGC 1399 and NGC 4472, 
respectively. We will refer to the latter by the globular cluster name RZ2109. The first has 
strong [OIII] and [Nil], the second only [OIII]. Both Ha and H)3 are very weak or undetected 
in both objects. We assume that the large line widths are due to Keplerian rotation around a 
compact object and derive expressions for maximum line luminosities. These idealized mod- 
els require central masses > 100 and > 30000 Mg for CXOJ03383 1.8-352604 and RZ2109, 
respectively. An independent, bootstrap argument for the total disc mass yields, for both sys- 
tems, Mdisc ^ 10^^ M© for a purely metallic disc (and two orders of magnitude larger for 
solar metallicities). If Roche-lobe overflow is implicated, viscous time-scales are > 300 yr. 
Standard disc theory then offers another limit on the central masses. Lobe radii for a ~ 1 Mq 
donor are > lO'^ cm. We therefore rule out Roche-lobe overflow of a white dwarf in both sys- 
tems. Red giants could fill the necessary lobes. Whether they are too metal-poor to produce 
the strong forbidden Unes without strong hydrogen emission is unclear 
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1 INTRODUCTION 

Ultraluminous X-ray sources are point sources that appear to ex- 
, ceed the Eddington luminosities of stellar-mass compact objects. 

■ They are therefore prime candidates for long-sought intermediate- 
I mass black holes (IMBHs). ULXs have been studied extensively 
. since their discovery in the early 1980s, and they continue to be 
' the subjects of intense research activity (Miller, Fabian, & Miller 
[ 2004; Winter, Mushotzky, & Reynolds 2006; Hui & Krolik 2008; 

■ King 2008, 2009). 

Zepf et al. (2007; 2008, hereafter Z08) and Irwin et al. 
(2010, hereafter 110) recently identified very broad and luminous 
[OIII] A5007 emission in two ULXs harbored by globular clus- 
ters in elliptical galaxies NGC 4472 and NGC 1399, respectively. 
The connections to globular clusters are robust, as in both 110 and 
Z08, we see large redshifts that agree very well with the redshifts 
of the clusters and their host galaxies. Line widths were 140 and 
1500 km s^' in the two systems, and line luminosities, assum- 
ing the sources radiate isotropically, were 10^^"" erg s-K Both 
Z08 and 110 found difficulty simultaneously explaining the line 
widths and luminosities with an accretion disc around a stellar 
black hole, the most natural explanation for a consistent source with 
Lj ~ lO^^'' erg s^ ' . The observed lines created a tension. Their large 
line widths suggested rapid rotation very close to the source, while 
the large luminosities suggested a large line-emitting region and an 
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origin in the outer disc. The lack of Ha and Hj8 presented further 
complications in the context of globular clusters. Z08 suggested 
a stellar black-hole wind is the explanation for RZ2109. 110 sug- 
gested for CXOJ03383 1.8-352604 tidal disruption of a white dwarf 
by an IMBH (Fabian, Pringle, & Rees 1975; Rosswog, Ramirez- 
Ruiz, & Hix 2009). 

This Letter is organized as follows. In section |2l we revisit 
and extend the arguments presented by Z08 and 110, first relax- 
ing the assumption that the line-emitting gas cannot be more dense 
than the critical density of the observed lines. We derive a simple 
expression for the maximum line luminosity in an accretion-disc 
geometry. Our expression is independent of temperature, density, 
and metallicity, and depends only on rotational velocity and cen- 
tral mass. While a super-critical density alleviates some of the ten- 
sion described above, we nonetheless confirm the Z08 and 110 con- 
clusions that stellar-mass accretion cannot be implicated in either 
ULX. In section [3] we present bootstrap arguments for the total 
mass in an accretion disc and obtain lower limits. In section|4l we 
explore the implications for Roche-lobe overflow. 



2 MAXIMUM LINE LUMINOSITIES 

Here we derive expressions for the maximum line luminosities of 
the [OIII] and [Nil] forbidden lines. We focus on the [OIII] lines 
because they are seen in both objects. Because both [Nil] and 
[OIII] are produced by ions of the carbon isoelectronic sequence. 
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Figure 1. [OIII] A5007 emissivity (red solid, erg cm^^ s^') and intensity 
(green dashed, erg cm^^ s"') emitted in an isotliermal, constant-density 
line-of-sight with solar metallicity and column density 1 cm"-. Each curve 
is a function of volume density. The critical density can be seen as a "knee" 



Figure 2. Theoretical ratio of the lines in the [Olll] doublet AA4959,5007 
as a function of column density. Observations in both Z08 and 110 
clearly suggest ratios toward the canonical optically-thin value. The vertical 
bar marks W(0+-)= 10" cm"^ (or T5007 = !)■ 



our formalism applies to both. The final expression is independent 
of temperature, density, and metallicity and depends only on rota- 
tional velocity and central mass. 

We assume the gas is photoionized and collisionally excited 
and begin by writing a simple two-level balance equation 



(1) 



where ni and 712 are the populations of the lower and upper lev- 
els, A21 is the spontaneous radiative transition probability (s"') 
from level 2 to level 1, and q\2 and ^21 are the coUisional exci- 
tation and de-excitation coefficients (cm-' s"' ) and are related by 
?12 = %i12\ exp{—hv/kT). The critical density Werit = ^2l/?2l is 
the density at which the collisional de-excitation rate («eg2l) equals 
the radiative decay rate. For <Si n„[i, the radiative decay domi- 
nates and the population of level 2 is 



"2 ■■ 



«c"l9l2 



(2) 



For rig 2> n^^i^, the usual Boltzmann equation of local thermo- 
dynamic equilibrium (LTE) applies. The local emissivity is e = 
nihViiAii, and it follows that 



neiiiqiihVii 



"e «^ "ciit 



"1 ^ exp(-/7V/fcr)/!V2lA2i ne > «crif 



(3) 



The critical density marks the transition from the regime where ra- 
diative cooling is proportional to ne njoj, to the LTE regime where 
cooling is proportional to n[g,^. The emitted energy per unit volume 
per unit time increases monotonically with increasing density. The 
two regimes are illustrated in Figure [T] for the [OIII] A5007 line. 
The critical density is represented by the knee near 10^ cm"'. 

Both 110 and Z08 considered the low-density case, and we 
have no need to reconsider that here. The remainder of this work 
will assume ne > /id-it. Authors who have discussed the high- 
density behavior of [OIII] emission include Nussbaumer & Storey 
(1981), Keenan & Aggarwal (1987), Kastner & Bhatia (1989), and 
Osterbrock & Ferland (2006). From an observational perspective, 
we note, for example, Andrea, Dreschel, & Starrfield (1994) found 
ne as high as 10^ cm"-' from [OIII] lines in classical novae. 

We generalize the high-density emissivity to multi-level con- 
figurations by using the 0+^ calculations of Nussbaumer & Storey 
(1981). Their calculations are valid for all practical densities and in- 
clude temperatures as high 40, 000 K, safely above the temperature 
of photoionized 0+^ gas considered in section[3]below. The upper 
level of the [OIII] A 5007 and [Nil] 16548 transitions is 'O2. The 
fractional population /('D2) in Table 5 of Nussbaumer & Storey 



does not exceed ~ 0.2. A comparable value applies to N^ 
emissivities take the maximum value 
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and our 



(4) 



where n^ is the density of the ionization stage. 

Next we posit a column of gas having local emissivity given 
by Equation|4] The intensity emitted through the column is 



edz, 



(5) 



where z represents position along the column (or above an annulus 
perpendicular to the plane of the accretion disc). To impose an up- 
per limit on the integration variable, we introduce a photon escape 
probability, /^sc = , where T is the (Napier) line-center optical 
depth. Equation[5]is then written 



(6) 



where we have used = n^Jdz. 

Equation [6] saturates at T ~ I, consistent with the rule-of- 
thumb that we can see into a cloud only up to optical depth unity 
(e.g., Rybicki & Lightman 1979). The upper limit corresponds to 
N{0+^) ^ 10^2 cm 2, assuming local line widths are thermal at 
^ K. Optical depth in the [OIII] lines causes their relative 
strengths to decrease from the canonical 3:1 ratio to approximate 
parity in the optically-thick limit. This behavior is shown in Fig- 
ure|2] The plotted ratio is visibly greater than unity in both Z08 and 
110. Our limit, therefore, has the consequence of preventing model 
line-ratios inconsistent with the observations. 

Note that Equation|6]does not depend upon density. The green 
dashed curve in Figure [T] illustrates this and is an important point 
in our analysis. Above the critical density, the intensity emitted 
through a fixed column density is independent of the volume den- 
sity of the gas. This means that, with no penalty on the total inten- 
sity, we can increase the density and, with all else fixed, squeeze 
the gas into a smaller volume, potentially allowing a stellar-mass 
black hole explanation. 

Finally, we obtain a maximum luminosity via an effective sur- 
face area. We define the inner radius, ro, of our line-emitting re- 
gion by assuming the line-emitting gas is in a Keplerian orbit about 
central mass M with velocity v (which we relate to observed line- 
widths below) so that 



ro = GM/v' 



(7) 



We assume the line-emitting region comprises an annulus with 
width comparable to rg. Line luminosities are written 



: 47trt,I 



on 



: 4k 



f GM 



(8) 
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Table 1. Observed and model luminosities (erg s ). The last column con- 
tains the values defined via equation [9] 

Quantity CXOJ03383 1.8-352604 RZ2109 



^5007 
^584 



few X lO^*" 1.4 X lO^' 2.4 x 10^3 
3 X lO^*" unobs. 2.8 x 10^2 

1.5-2.3x1039 4x103'' 



We evaluate constants, normalize variables to convenient values, 
and obtain maximum luminosity 



?0m2„ -4 
■^x'^l I'lOO • 



(9) 



where Mi is the central mass in units Mq, vioo is the rotational 
velocity in units 100 km s^' , and all constants have been absorbed 
into fiducial values . Table [T] contains the fiducial luminosities 
for [OIII] 15007 and [Nil] 16584. Atomic data used in this section 
are from Wiese, Fuhr, & Deters (1996), as obtained from the NIST 
Atomic Spectra Database (Ralchenko et al. 2008). 

It is important to emphasize that equations [8] and |9] are inde- 
pendent of temperature, volume density, and metallicity. This is by 
design. However, temperature and metallicity do inform the prac- 
ticality of the equations. We will discuss the former below. Re- 
garding the latter, at fixed n^, the disc height required to reach a 
certain column density (and therefore luminosity) is inversely pro- 
portional to metallicity. A purely metallic plasma allows the small- 
est disc height for the maximum-luminosity configuration of our 
model. Solar metallicity gas requires a disc height more than three 
orders of magnitude greater, reaching ~ 10 Thomson depths. Scat- 
tering would yield emergent line widths --^ 1000 km s^^. Such 
large widths would contradict our rotational broadening assump- 
tion in both ULXs and are ruled out completely by observation in 
CXOJ03383 1. 8-352604. If the width scales as the number of scat- 
terings (usually the greater of r and T^), this problem would be 
eliminated by decreasing the column by some factor. The mini- 
mum central mass found below would increase by the square of 
that factor. 

We plot in Figure [3] maximum line luminosities as a function 
of the velocity of the gas at tq. Several values of the compact ob- 
ject mass are considered. The middle two lines represent recent 
theoretical (80 Mp^) and observed (30 M,;;) upper limits to the mass 
of a stellar black hole (Belczynski et al. 2010). Results for 110 and 
Z08 are indicated by the squares, where we have assumed veloci- 
ties equal to l/(2\/ln2) times the measured FWHM. This ignores 
the effect of disc inclination on the observed line widths. Account- 
ing for inclination would mean a greater rotational velocity and a 
smaller annulus radius (for a given central mass). This would de- 
crease the surface area of the line-emitting region and the lumi- 
nosity of the line. The arrow on each square points to the velocity 
appropriate for the (arbitrarily chosen) inclination angle 45 ° . 

The point corresponding to CXOJ033831. 8-352604 in the up- 
per panel of Figure[3]is beneath the lines corresponding to the theo- 
retical and observed stellar-mass black hole upper limits. However, 
the lower panel of Figure [3] makes a more compelling case. The 
observed L6584 is very slightly less than (at edge-on inclination) 
the maximum predicted luminosity corresponding to the maximum 
theoretical mass of a stellar black-hole (Belczynski et al. 2010). 
Note that these predictions assume optimal cooling efficiency. 110 
found Tg ^ 13,000 K in the [Nil] emitting region using the lack of 
15755. Using this temperature would lower the [Nil] predictions 
by a factor ~ 2 and require M > 100 M0 . [OIII] 14363 was neither 




Figure 3. Maximum line luminosities (upper: [OIII] A5007; lower: 
[Nil] A6584) versus rotational velocity of the line-emitting region for sev- 
eral values of the central mass (in units M, . ). Squares indicate results from 
110 and Z08. Arrows point from minimum (edge-on) velocity to the velocity 
appropriate for 45° inclination. [Nil] was not detected in RZ2109. 



presented nor mentioned in 1 10. Measurement of this line could 
place a strong constraint on the temperature of the [OIII] region. 

The upper panel of Figure[3]clearly demonstrates that the Z08 
observations of RZ2109 cannot be explained as Keplerian rotation 
around a stellar-mass object. The required mass is at least 6000 Mq. 
[OIII] 14363 is weakly detected in RZ2109. A rough estimate is 
that 15007 is at least 10 times greater. The Nussbaumer & Storey 
(1981) tabulations would then give T < 15,000 K at the 15007 
critical density. That upper limit quickly decreases to T < 7, 500 K 
for "e }t 10^ cm^-'. These considerations depress the maximum 
luminosity by at least a factor of 20, so that the minimum central 
mass is at least 30,000 Mq. 

We rule out stellar black-hole Keplerian accretors in both 
sources. Sub-Keplerian motions have orbiting radii that are smaller 
for a given central mass than Keplerian, decreasing the surface area 
of the line-emitting region and exacerbating the problems discussed 
above. We can therefore rule out sub-Keplerian, stellar black hole 
accretion more readily than we can rule out Keplerian. 

We stress that the two panels of Figure[3]do not depend upon 
the relative abundances of oxygen and nitrogen. The predicted 
[Nil] luminosities are noticeably less than the [OIII] luminosities at 
the same velocity and central mass. This is because the [Nil] lines 
are less efficient coolants than their [OIII] counterparts. 



3 TOTAL DISC MASS 

Here we present a bootstrap argument for the total disc masses in 
an accretion scenario. The mass required to emit L5007 (or L^5S4) is 



MCD2): 



Zm„L 



hvA 



(10) 



where Zirip is the mass of a single atom with mass number Z. This 
is independent of the geometrical model and the gas density (even 
below n^nJ, and accounts only for the population in the upper level 
of the transition. Observed values of L^^ are listed in Table [T] The 
[OIII] lines require masses 3 x 10^^ and ~ 2 x 10^^ Mq for 
CXOI03383 1. 8-352604 and RZ2I09, respectively. The [Nil] lines 
in CXOI03383 1.8-352604 require 4 ~ 10"*' Mq. 

We can bootstrap our way to a total disc mass via 
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Figure 4. Ionization fractions of oxygen as a function of tlie temperature of 
an intense blackbody. 
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Figure 5. Ionization fractions of oxygen (upper panel) and nitrogen (lower 
panel) as a function of the electron temperature of gas ionized by a power- 
law continuum with spectral index F = 2.5. 



M(l£)2) _ Pz IfionfCD2)rdr 



Mdisc 



/ rdr 



(11) 



where /(ID2) was discussed in section|2]and is the C-like ion- 
ization fraction, which we estimate using the well-known plasma 
simulation code Cloudy (version COS, last described by Ferland et 
al. 1998). Figure |4] plots the results for oxygen as a function of 
the temperature of a bright blackbody photoionization source. The 
fraction in 0+~ peaks around 20,000 K. The maximum equivalent 
width of A5007 with N{0+^)= 10^^ cm"^ against a blackbody of 
the same area and Tgg = 20,000 K is about lA. We estimate from 
spectra in 110 and Z08 observed equivalent widths of 5 and 15 A. 
Blackbody sources with weak enough continuum at 5000A are 
too weak to ionize 0+. We require a harder continuum. 

We instead consider, again using Cloudy, a power-law contin- 
uum with spectral index T = 2.5 (following X-ray continuum fits 
in 110). Figure |5] plots both oxygen and nitrogen ionization frac- 
tions as a function of electron temperature. The ionization parame- 
ter ranges from — 9 ^ logC ^ 0. 

The fractions /jo,, and /('£)2) depend on temperature. We re- 
late them to the variable of integration r via the standard Shakura- 
Sunyaev prescription (T «: r"^/"*). The ratio of integrals in equa- 
tion [TT] is illustrated by the dotted blue curves in Figure |6] The 
result depends on the upper limit of integration. The peak is ~ 4% 
for both oxygen and nitrogen. The curves decline as r^^ at large 
radii. 

An important conclusion can be derived from Figures [5] and 
(6] The regions emitting [OIII] and [Nil] are mostly non-cospatial. 
The details depend on the ionizing source and the geometry. Since 




ratjius (arbitrary urrits) 



Figure 6. Level and ionization stage fractions of oxygen (upper panel) and 
nitrogen (lower panel) versus normalized radius. Solid red lines are frac- 
tions in O^^ and N^, dashed green lines are /('D2), and dotted blue lines 
are ratios of integrals given in the right-hand side of eauation fTTI 



the second ionization energy of nitrogen is only 85% of the second 
ionization energy of oxygen, the basic result is fairly robust. N+ 
always peaks at a temperature significantly less than the tempera- 
ture of the O"'"^ peak. If the standard a-disc is applicable, we have 
T oc r^-'/'*. Keplerian rotation gives v o= r^'/", so v T^^^. A se- 
ries of tests suggests the [OIII] and [Nil] line widths should differ 
by a factor > 1.4. The widths are uncertain but appear compara- 
ble in the 110 observations. This argues against Keplerian rotation 
about any central mass in CXOJ03383 1.8-352604. A very steep 
temperature gradient or a thermal instability between the [OIII] and 
[Nil] regions could mitigate this problem. High-resolution spec- 
troscopy would be useful in efforts to address such questions. 

Finally, we require the density of oxygen relative to the total 
gas density. Assuming metallicity is independent of radius, we can 
write 



Pz 
P 



Anz 



(12) 



I. Anz 

where A is atomic mass number. For fixed relative metal abun- 
dances, the fraction of mass in oxygen varies from ~ 0.005 with 
solar abundances (Grevesse & Sauval 1998) to ~ 0.4 for all metals. 

Combining the above results, we find Mdisc > 10""^ Mq. This 
corresponds to the purely metallic case. If we assume roughly so- 
lar abundances, Mjisc > 10^^ Mq. An estimate for the upper limit 
is not readily apparent in either system because we have no way 
of estimating the outer radius and, in particular, because the line- 
emitting gas could be only a thin outermost layer of the disc. 



4 ROCHE-LOBE OVERFLOW 

In steady Roche-lobe overflow accretion, the disc mass must be 
replenished on the viscous time-scale 



M 



(13) 



where the mass-accretion rate M = L^/rjc^ and T] is an efficiency 
typically taken to be 0.1. We can also consider the viscous time- 
scale in a standard a-disc. We manipulate the familiar expression 
(Frank, King, & Raine 2002), substituting azimuthal velocity for 
the radius via equation [T] and obtain 
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W ^ 4.2a-4/5M-//iV,3/2v,oo-^/2 yr, (14) 



where Mig is the mass accretion rate in units lO'^g s For the 
X-ray luminosities given in Table [T] we find M ~ 2 x lO'^g s^' 
(3 X 10"^ Mq yr"^). 

If we assume the "edge-on" velocities discussed in relation to 
Figure[3]and a = 0.1, equations 1131 and [74l require masses 70 and 
4000 Mq for CXOJ03383 1.8-352604 and RZ2109, respectively. 
These numbers are surprisingly similar to the minimum masses 
found in Section |2] The minimum masses are --^ 40 times larger 
if the disc has Z ~ O.3Z0. The viscous time-scale /vise ^ 300 yr. 

Temperatures derived with the a-disc prescription, standard 
assumptions, and the above limits are too cold by a factor ~ 10. 
Ml , M, and a all enter as fairly weak powers and offer no obvi- 
ous solution. Irradiation of a surface layer (e.g., Dubus et al. 1999) 
might yield the necessary temperatures. 

Can a globular cluster form and sustain the massive accretion 
discs implied by the observed [OIII] luminosities? Combining the 
minimum masses obtained in Section [2] with the rotational veloc- 
ities derived above, we obtain ro ~ 2 X 10'"* and 5 x lO'"* cm, for 
CXOJ03383 1.8-352604 and RZ2109, respectively. We take this as 
the distance from the center of the compact object to the Lagrange 
point Li . We then assume a 1 M,;; donor star (giving mass ratios 
q = 0.01 and 1.7 x 10^"*) and employ approximations from Eggle- 
ton (1983) and Frank, King & Raine (2002), and we obtain donor 
Roche lobe radii ri > lO'-' and 5 x lO'^ cm, respectively. These 
radii are comfortably within the size limits of the red giants that 
are common in globular clusters, but they represent lower limits. 
Larger central object masses would correspond to larger overflow 
radii (because the dependence of r^/ro on M via q is not enough 
to overcome ro = M). White dwarfs have been mentioned because 
of the apparent high-metallicity in both systems. With prototypical 

= 0.6 M0 and r„.j = 10^ cm, we can strongly exclude white 
dwarfs as Roche-lobe overflow donors in either system. 

Does the envelope of a red giant contain enough oxygen to 
produce the observed line luminosities? Carretta & Gratton (1997) 
considered the metallicities of red giants in 24 galactic globular 
clusters. They found sub-solar metallicities for every star in their 
study. The majority were sub-solar by at least an order of magni- 
tude. The most metal-rich had Z ~ 0.3Zp;. 

The [0III]/Hj8 and [NII]/Ha values 110 observed in 
CXOJ03383 1.8-352604 (> 5 and > 7, respectively) and the 
[OIII]/Hj3 value Z08 observed in RZ2109 (~ 30) are not neces- 
sarily indicative of high-metallicity gas. The photoionized gas near 
the center of the giant HII region 30 Doradus has Z ~ O.3Z0 and 
[OIII]/H/3 ~ 6 (Pellegrini et al., submitted). Elevated ratios are of- 
ten thought to be signs of radiative shock heating (e.g., Dopita & 
Sutherland 1995), as was noted by Z08. Detailed attempts to dis- 
tinguish between enhanced abundances and non-equilibrium heat- 
ing processes are probably unwarranted without better observa- 
tional constraints. A further complication pertains to whether the 
line-emitting gas is radiation- or matter-bound (McCall, Rybski, & 
Shields 1985), although the high observed [NII]/Ha would seem to 
make the latter less likely in CXOJ033831. 8-352604. 

If Z08 are correct in their black-hole wind explanation of 
RZ2109, we offer without comment the suggestion that the ob- 
served very broad line-widths may be the low-resolution appear- 
ance of narrow emission lines atop broader ones as in a two-wind 
structure (e.g., Fernandes 1999). 
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